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Abstract 

We present an experimental and theoretical investigation of the Generalized Ein- 
stein Relation (GER), a particular form of a fluctuation-dissipation relation, in an 
out-of-equilibrium visco-elastic fluid. Micrometer beads, used as thermometers, are 
immersed in an aging colloidal glass to provide both fluctuation and dissipation 
measurements. The deviations from the Generalized Einstein Relation are derived 
as a function of frequency and aging time. The observed deviations from GER are 
interpreted as directly related to the change in the glass relaxation times with aging 
time. In our scenario, deviations are observed in the regime where the observation 
time scale is of the order of a characteristic relaxation time of the glass. 
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1 Motivation 



During the past years, slow relaxation systems have received considerable at- 
tention in the sense that they constitute a challenge for new developments in 
nonequilibrium statistical mechanics. Such systems, which are rather common 
in nature, include structural glasses, polymers, colloids and granular matter. 
They share some characteristics called glassy behavior, in particular a drastic 
slowing down of relaxation processes when some control parameters are varied. 
As the characteristic relaxation times becomes longer than, or comparable to, 
the observation timescale, the system is said to age. The physical properties 
of the material depend on the waiting time, which is the time elapsed since 
preparation in the nonequilibrium state. 

Recently, there have been serious attempts to approach and understand ag- 
ing phenomena on a theoretical level. Many efforts have been devoted to ap- 
ply statistical physics description - in particular, fluctuation-dissipation rela- 
tions (FDR) - to out-of-equilibrium systems. At thermodynamic equilibrium, 
fluctuation-dissipation relations relate the response functions of the system to 
its autocorrelation functions, involving a single thermodynamic parameter, the 
equilibrium temperature. However, FDR are not expected to hold in out-of- 
equilibrium systems. The idea that a timescale dependent nonequilibrium or 
effective temperature could describe the slow relaxation modes has received a 
lot of attention. This effective temperature, defined from an extension of FDR, 
is different from the bath temperature and has been shown to display many of 
the properties of a thermodynamic temperature yj; Deviations from FDR 
have been observed in many numerical simulations (3; 4; 5; 6; 7; sj; oj; 10). To 



date, there are still few experiments where FDR are studied in aging mate- 
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rials, such as structural glasses dill), c olloidal glasses yj 



spin glasses (18), polymer glasses 
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and granular media (120 
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In most 



cases, they show deviations from FDR as a function of waiting time and fre- 
quency. Unfortunately, no global understanding can yet be extracted from 
these experimental findings. 



Here we present an experimental and theoretical investigation of the Gener- 
alized Einstein Relation (GER), a particular form of a fluctuation-dissipation 
relation, in an out-of-equilibrium visco-elastic fluid. The deviations from the 
GER are measured as a function of waiting time and frequency u, leading 
to an effective temperature Teff(cij), parametrized by t^. This is achieved by si- 
multaneously measuring the response function to an external force - using an 
optical tweezer - and the position fluctuations of micrometric beads embedded 
in the glass. 



The paper is organized as follows. Sec. [2] is devoted to the theoretical investi- 
gation of the Generalized Einstein Relation in an out-of-equilibrium system, 
leading to the definition of an effective temperature. Sec. [3] presents the ex- 
perimental procedures used to test such a GER in an aging colloidal glass: 
Brownian motion measurements and response to an oscillatory force - using 
optical tweezer. In Sec. HI the effective temperature is derived as a function of 
aging time and frequency. In Sec. [5l we suggest a microscopic interpretation 
of our findings. We conclude the paper with comparisons with other tests of 
FDR in the same glass and with open questions. 
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2 Fluctuation-dissipation relation in an out-of-equilibrium envi- 
ronment 



2.1 Equilibrium medium case 



Let us first consider a diffusing particle of mass m evolving in a stationary 
medium. Its motion can be described by a generalized Langevin equation, in 
which v{t) = dx/dt is the particle velocity, F{t) the random force acting on 
the particle, and 7(t) a delayed friction kernel that takes into account the 
visco-elastic properties of the medium : 

m^ + m 1 ^{t-t')v{t')dt' = Fit) (1) 

— oo 

Note for further purpose that, in the usual experimental frequency range, the 
inertial term in Eq.([T]) is quite negligible. The frequency- dependent particle 
mobility thus reduces to jl{uj) = [m^{uj)]^^, where 7(1^) = 7(^)e*'^* is 
the Fourier transform of the delayed kernel 7(t). 



2.1.1 The Generalized Einstein Relation 

If the surrounding stationary medium is in thermal equilibrium at temperature 
T, one can write a Kubo formula expressing the frequency- dependent particle 



mobility fi{w) in terms of the velocity correlation function (l22l ) : 



fiiiu) = ^ / {vit)vit = 0))e-*dt (2) 







(the symbol (. . .) denotes the average over an ensemble of realizations). For- 
mula ([2]) defines fi{uj) as an analytic function of uj in the upper complex half 
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plane. 

From Eq. ([2]), one can derive a formula linking the drift and diffusion properties 
of the particle. The mean-squared displacement of the diffusing particle, as 
defined by : 

{Ax^{t)) = {[x{t) -x{t = 0)f), t>0 (3) 

can be deduced from the velocity correlation function via a double integration 
over time: 

t ti 

{Ax\t)) = '^jdhj dt2{v{h)v{t2)) (4) 



Introducing the Laplace transformed quantities v{s) = v{t)e~^^dt and 
(Ax^(s)) = {Ax'^{t))e~^^dt, one gets, by Laplace transforming Eq. ([3]): 

s^A^^is)) = 2{v{s)vit = 0)) (5) 

The quantity {v{s)v{t = 0)) can be obtained by setting uj = is in the Kubo 
formula ([2]) (which is allowed for positive s, since the function il{ui) is analytic 
in the upper complex half plane). One gets : 

{v{s)v{t = 0)) = kT fx{s) (6) 

with fi{s) = jl{uj = is). Eq. ([5]) then takes the form of a relation linking the 
Laplace transforms of the particle mean-squared displacement and mobility : 

s^A^^s)) = 2kT fi{s) (7) 
Eq. ([7j) is known as the Generalized Einstein Relation (GER). 
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2.1.2 The fluctuation- dissipation relation 

It follows from Eq. ([2]) that the velocity spectral density is related to the 
dissipative part of the mobility by a fluctuation-dissipation relation (namely, 
the celebrated Einstein relation): 



The validity of Eq. ([8]) is restricted to uj real. 
2.1.3 Discussion 

The Einstein relation ([H]) contains the same information as the GER ([7]), since 
jl{uj) can be deduced from 3?e/i(c<j) with the help of the usual Kramers-Kronig 
relations valid for real u. Thus, Eqs. ([7]) and ([8]) constitute fully equivalent 
formulas, which both involve the thermodynamic bath temperature T. 

2.2 Out- of -equilibrium environment 

The general situation of a particle diffusing in an out-of-equilibrium environ- 
ment is much more difficult to describe. As well-known, in an aging medium, 
no well defined thermodynamical temperature does exist, so that Eqs. ([7]) and 
dH]) are no longer expected to be valid. 

We address the question whether the study of the diffusion and drift of the 
probe particle is likely to provide information about the out-of-equilibrium 
properties of its surrounding medium. 



oo 
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2.2.1 Modified Generalized Einstein Relation 

Out of equilibrium, one can write a modified Kubo formula as 



/iM = -^ /MtMt = 0))e-*rft (9) 







Formula ([9]) defines the product /i(u;)0(u;) as an analytic function of uj in the 
upper complex half plane. Since this analyticity property holds for /i(a;), it 
also holds for Qioj). 

Accordingly, Eq. ([6]) has to be replaced by : 

= 0)) = A;e(s)/i(s) (10) 



Thus, the relation linking the Laplace transforms of the particle mean-squared 
displacement and mobility writes : 

s^{Ki^{s)) = 2kQ{s)fi{s) (11) 

Eq. flTT]) states an out-of-equilibrium GER. It involves a s-dependent GER 
ratio 0(s), parametrized by the age of the system tw. 

Experimentally, it is always possible to measure independently both quantities 
- (Ax2(s)) and /t(s) - during aging. For a particle evolving in a bath at 
equilibrium, 0(s) = T would be independent on s, as shown by Eq. ([7]). In an 
out-of-equilibrium environment where G)(s) is a priori expected to be different 
from T, the mobility of the probe particle cannot be deduced from the mean- 
squared displacement measurement, and vice-versa. As a consequence, the 
visco-elastic properties of the out-of-equilibrium medium cannot be deduced 
from the single measurement of the probe thermal fluctuations and passive 
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microrheology cannot be achieved [21]. 



2.2.2 Modified fluctuation-dissipation relation 



Since, in an out-of-equilibrium environment, even stationary, the Einstein re- 
lation ([8]) is not a priori satisfied, it has been proposed yj; |2|) to rewrite it in 
a modified way with the help of a frequency-dependent effective temperature. 
Such a quantity, denoted as Tce{ijj) and parametrized by the age of the system 
t^, has been defined, for real tu, via an extension of the Einstein relation ([8]). 



One writes a modified Einstein relation as (124 



m ■ 



(12) 



It has been argued in ([l|; 12) that the effective temperature Tes{uj) defined 
in such a way would have a possible thermodynamic meaning, in the sense 
that it plays the same role as the thermodynamic temperature in systems 
at equilibrium (namely, it controls the direction of heat fiow and acts as a 
criterion for thermalisation). 



2.2.3 Discussion 

In the case of an out-of-equilibrium medium, one has at hand the modified 
GER ffTTl) and the modified FDR f|T2|) . The question thus arises of the link 
between these two descriptions, namely of the relation between the GER ratio 
6(s) and the effective temperature Teg (a;). 

The modified Kubo formula ([9]) defines a complex frequency-dependent func- 
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tion Q{uj), in terms of which the velocity spectral density writes: 



C^^icu) = 2me[e{uj)il{uj)] (13) 



Comparing Eqs. f[T21) and f[T5]l . one gets the following relationship between 
Q{lj) and Tes{uj) : 

T,s{uj) ^ej2{uj) = ^e[j2{uj)e{uj)] (14) 



Eq. f|T4l) displays the fact that the effective temperature Tcs{uj) involved in the 
modified Einstein relation ( |T2l) can be deduced from the GER ratio introduced 
in the modified GER (fTTl) (and vice-versa) . This point has been developed in 



details in (l25l ) 



Thus, as it can be seen from Eq. (|TT|) . independent measurements of the parti- 
cle mean-squared displacement and mobility in an aging medium give access, 
once (Ax^(s)) and fi{s) = jlioj = is) are determined, to the GER ratio 6(s). 
Eq. (HM then allows to derive an effective temperature Teg (a;) from the ex- 
perimentally measurable GER ratio Q{uj) = G)(s = —iu). 

In the out-of-equilibrium case, the modified GER ffTTl) and the modified FDR 
( |T2i) constitute fully equivalent formulas. However, Eq. ( fTTl) involves a GER 
ratio 0(s), naturally expressed in terms of the Laplace variable s, while Eq. 
(fT2l) involves an effective temperature Teff(u;), naturally expressed in terms of 
the real frequency u. 
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3 Experimental procedures 



3. 1 Samples preparation 



The experiments were performed on Laponite RD, a syntlietic clay manufac- 
tured by Laporte Industry. The particles of Laponite are colloidal disks of 25 
nm diameter and 1 nm thickness, with a negative surface charge on both faces 



(l26l) ■ The clay powder was mixed in ultra-pure water, and the pH value of 



the suspensions 



ixed to pH = 10 by addition of NaOH, providing chemically 



stable particles (l27l). The suspension was stirred vigorously during 15 minutes 
and then filtered through a Millipore Millex - AA 0.8 /im filter unit. This 
procedure allows us to prepare a reproducible initial liquid state. The aging 
time = of the suspension is defined as the moment it passes through the 
filter. 



hese aqueous suspensions form glasses for low volume fraction in particles 



((281). Starting from a viscous "liquid" state right after preparation, the sus- 
pension becomes more and more visco-elastic with time. Since the physical 
properties of the suspension depend on the time elapsed since preparation, 
the sample is said to age. Aging can be seen through the chan ge i n both the 



30l ). Laponite 



visco-elastic properties and of the colloidal disks diffusion (129 : 
suspensions age on timescales that depend on the particles' concentration. We 
are thus able to control the aging timescales of the glass by adjusting this con- 
centration. With a volume fraction of 2.3% wt, the glass evolves over several 
hours, slowly enough to allow us quasi-simultaneous measurements of the fluc- 
tuation and dissipation properties, within a few minutes, without significative 
aging of the sample. These two successive measurements are thus considered 
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to be performed at the same waiting time tyj. 

The experiments were carried out in a square chamber - 20 x 20 mm^ - made 
of a microscope plate and a coverslip separated by a thin spacer (0.1 mm 
thickness). The beads are suspended in the glass right after its preparation. 
The chamber is then filled with the suspension, sealed with vacuum grease and 
mounted on a piezoelectric stage on the optical microscope plate. The probes 
are latex and silica beads, in very low concentration (respectively 10~^% and 
4.10~^% in volume). Latex beads (1.0±0.1//m in diameter, Polysciences, Inc.), 
were preferentially used for fluctuation measurements : since they do not sed- 
iment during the experiment, their random motion is not perturbed by the 
chamber walls. Sihca beads (2.1 ± 0.1/im in diameter. Bangs Lab Inc.) were 
used for dissipation measurements, because they are more efficiently trapped 
by the optical tweezers. By assuming that the particle mobility scales as the 
inverse particle size according to fl{uj) = l/Q%Rfi{ijj) (Stokes assumption in a 
visco-elastic fluid) , the comparison between the results of the fluctuation and 
dissipation measurements, once rescaled to the same diameter, is meaningful. 
This relation is exact in a purely viscous fluid at low Reynolds number; here 
we assume its validity in the visco-elastic fluid. 

In these experiments, the nature of the bead material may play a role through 
the particle-sample interactions, dominated by electrostatics. Actually, the 
suspensions of Laponite are prepared at pH = 10 by addition of NaOH. At 
this ionic concentration, the Debye-Hiickel screening length is of the order of 30 
nm, which is comparable to the Laponite particle diameter, and smaller than 
the interparticle distance. In such conditions, the electrostatic interactions 
between the micro beads (actually both sihca and latex beads are charged) 
and the Laponite particles are screened on a distance much smaller than the 
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bead diameter. As a consequence, the range of interactions between beads 
and particles does not exceed 30 nm and the system is not perturbed on large 
scales by the introduction of the beads. 

3.2 Thermal fluctuations of the beads 

At a given aging time i^, we record the fluctuating motion of 1 /im latex beads 
during 8 s, with a fast CCD camera sampling at 250 Hz (Fastcam-PHOTRON 
LTD). A digital image analysis allows us to track the bead positions x{t) and 
y{t) close to the focus plane of the microscope objective. For each bead, we 
calculate the time-averaged mean-squared displacement {Ar'^{t))t' = {[x{t' + 
t) - x{t')]'^ + [y{t' + t)- y{t')f)t' = 2{Ax'^ {t))t'. To preserve a good statistics, 
we keep the data of {Ar'^(t))t' in the range 0.004 < t < 1 s. The resolution 
on the bead position, determined by sub-pixel accuracy in the image analysis 
detection, is 10 nm. The glass remains in a quasi-stationary state during the 
recording, which takes a short time compared to the aging timescale. The 
quantity {Ar^(t))t', averaged over several beads and realizations, can thus be 
identified to the ensemble-averaged mean-squared displacement. 

3.3 Dissipative response to an oscillatory force 

We describe now the measurement of the mobility at a given frequency 

cu, for various aging times tw This measurement is performed immediately 
after the Brownian motion recording, at the same aging time tw Since the 
aging Laponite suspension is a visco-elastic fluid, the bead mobility jl{u!) — 
|/2(a;)|e*'^^'^^ is a complex number. We thus need to measure the phase and 
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modulus of the tracer mobility. We use an optical tweezer to trap a 2.1 /im 
silica bead immersed in the glass. Trapping is achieved by focusing a powerful 
infrared laser beam (Nd YAG, Spectra-Physics, Pmna: = 600 mW) through a 
microscope objective of large numerical aperture (jsil). The trapping force F 
on a small dielectric object like a silica bead is proportional to the intensity 
gradient in the focusing region. It depends on the distance x of the center of 
the trapped object to the center of the trap, according to F{x) = — A;x(l+e|x|). 
The corrective factor e is introduced to take into account the non-harmonicity 
of the trapping potential. The trap stiffness k and the factor e are determined 
from an independent calibration. The calibration procedure was described in 



details in previous publications using optical tweezers set-up (l34j ). Once the 
bead is trapped, we make the experimental chamber oscillate by monitoring 
the displacement Xpexp{—iut) of a piezoelectric stage. Neglecting in a first 
step the visco-elastic fluid inertia, the fluid displacement xj in the bead vicin- 
ity is taken equal to the piezolectric stage displacement Xp. The validity of this 
assumption is discussed in the Appendix. Due to the relative bead/fluid mo- 
tion, the visco-elastic fluid exerts a sinusoidal force F' exp{—iujt) on the bead. 
We record with the fast camera the bead movement, and measure by con- 
ventional image analysis its displacement xexp{—iujt) from the trap center. 
Notice that x is a complex number which includes a phase shift due to the fluid 
visco-elasticity. At a given frequency u, the force amplitude F'{uj) is given by 
F'{lj) = v{ijj) / where v{ijj) = iuj{xf—x) = iuj{xp—x) is the relative cham- 
ber wall / bead velocity, and fi{uj) the Fourier transform of the bead mobility. 
In our range of experimental frequencies (0.5 < / < 10 Hz), the bead inertia 
is negligible, so that we can simply use the relation F + F' = to calculate 
|/i(co')| and f{uj). Notice that the motion of the piezoelectric stage is numeri- 
cally controlled by a sequence of successive sinusoidal signals at five different 
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frequencies {0.5, 1, 2, 5, 10} Hz. The same program synchronously generates a 
sequence of pulses to trigger the image acquisition, so that the phase shift be- 
tween the force and the bead movement can be accurately measured. For both 
the passive and the active microrheology experiments, we used a conventional 
image processing software (IMAQ vision builder from National Instruments) 
to determine the bead position from the numerical recording of the bead im- 
ages. The accuracy on the bead position was about 10 nm, for typical motions 
comparable to the bead radius, in the micrometer range. 

4 Results 

4-1 Rescaling of the aging dynamics 

Figure [T](a) shows the increase in the viscosity modulus with waiting time 
for a set of six different realizations. The complex mobility of a Laponite 
suspension at 2.3% wt was measured by applying an oscillatory force on a 
silica bead in the linear regime. The complex viscosity modulus |^(c(;)| was 
estimated from the complex mobility modulus 1/2(^^)1 assuming that the Stokes 
relation fl^uj) = ^^^^-^^^ remains valid in the visco-elastic fluid. Starting from 
a value close to the water viscosity 10^'^ Pa.s, the glass viscosity modulus is 
shown to increase by 3 orders of magnitude over about 200 — 400 minutes. 
These variations in aging timescales are not so surprising because in such 
suspensions, the aging dynamics is known to drastically depend on the particle 
concentration. In a suspension at 2.5% wt, a three orders of magnitude increase 
of the complex viscosity modulus - measured with a conventional rheometer - 
is reached after 100 minutes, while in a 3.5% wt suspension, the same increase 
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is reached after 10 minutes (3(f). Besides, slight ionic concentration variations 
in the suspension induce important changes in particle interactions. These 
variations may arise either from pH differences or ions salting-out from the 
chamber walls. 



Here we suggest a time rescaling in order to compare the different realizations. 
We assume that the differences in the aging dynamic shown in Fig. 1(a) can 
be compensated by a linear stretching in aging time. This is supported by 
a recent paper by Joshi where the glass transition of a Laponite suspension 



is investigated (1321 ). The characteristic dimensionless relaxation time t/tq of 
the suspension before the system enters the full aging regime {t^ < t^*) is 
shown to only depend on the dimensionless age t^/t^* . The transition time 
tw* depends on the volume fraction $, according to t*^ ~ $"2/^. Given a mas- 
sic fraction of 2.3 %, and the aging timescales involved, our experiments were 



performed before the full aging regime takes place (1291 ). By assuming that i) 
the suspension viscosity only depends on the characteristic relaxation time 



t/tq and thus on t^jtw* (1331 ) . and ii) that the differences between realizations 
are essentially due to slight volume fraction variations, it seems realistic to 
balance the differences by a linear stretching in aging time. In Fig. [I](b), the 
aging time was linearly stretched to make the viscosities roughly coincide for 
all experiments, choosing a particular realization as a reference. Nevertheless, 
it must be clear that the same conclusions on T^q are obtained on each sep- 
arate sample. With this rescaling, we are now able to average over different 
realizations. From now on, the aging times of the different realizations are 
rescaled by using this rescaling. 
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Fig. 1. (a): Viscosity modulus of the glass as a function of aging time, for six different 
realizations. The viscosity modulus was measured by applying an oscillatory force 
at a frequency / = 1 Hz on a silica bead in the linear regime; (b): Rescaling of the 
complex viscosity modulus. The aging time is linearly stretched t?^ — > = atw 
to make the viscosities roughly coincide for all experiments, choosing a particular 
realization as a reference. 

4-2 Brownian motion of heads in the glass 

Figure [2] shows the mean-squared displacement of latex beads immersed in the 
colloidal glass, as a function of time t, for various aging times tw At t^, = 0, we 
observe a nearly diffusive behavior of the tracer beads, characterized by a linear 
dependency of the mean-squared displacement with time. Upon increasing on 
tw, the tracer motion becomes sub-diffusive. 
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Fig. 2. Mean-squared displacement of 1 /im latex beads immersed in the 
glass, as a function of time. The curves correspond to different aging times 
tm = 15, 64, 100, 148, 169 and 186 minutes from top to bottom, measured for one 
realization. The fluctuating motion is purely diffusive at short t^,, and becomes 
sub-diffusive as the glass ages. 
4-3 Mobility from oscillatory force 

Figure [3] shows the frequency and aging time dependence of the complex mo- 
bility modulus and phase (p{uj) measured when applying an oscillatory 
external force at successive frequencies {0.5, 1, 2, 5, 10} Hz. The mobility mod- 
ulus |/i(ci;)| of the bead was found to decrease with aging time t^ as shown in 
Fig. [3](a). This corresponds to the increase in the glass visco-elastic modulus. 
Fig.[3](b) represents the phase shift (p{uj) as a function of t^. Starting from zero, 
the phase ip decreases with t^. This corresponds to the shift between a purely 
viscous behaviour = 0) to a visco-elastic one (— 7r/2 < (/? < 0). In a first- 
order approximation, the phase remains independent of the frequency. Figure 
[2(c) shows the same data, plotted as a function of the frequency u, for 
different aging times t^. The mobility modulus |/i(co')| is well fitted by a power 
law = HoiuJ /ojqY in the experimental frequency range corresponding to 
one and half decade. Here, ujq is an arbitrary reference frequency set to = 1 
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rad.s"^ for convenience. Starting from about zero at low aging times, the ex- 
ponent (3 increases with t^. Since |/i(c(;)| exhibits a power-law behavior, the 
dependence of ^{t^) must be related to [3 hy ipity^) = — /5(t^)7r/2. This ensures 
that jl{s) is a real number. In Fig. [3](b), we plot the quantity — /5(t^)7r/2 (full 
circles), which is in fair agreement with the experimental data f{tw)- 

4-4 Effective Temperature derived from oscillatory measurements 

The GER ratio G)(s) is now derived from the oscillatory and Brownian mo- 
tion measurements. The selected range s G [3; 60] s~^ corresponds to the 
experimental frequencies range / G [0.5; 10] Hz. From the oscillatory mea- 
surements, the variations of fi{uj) are well represented by the analytical form 
/i(co') = fio^u/uo)'^ exp {—iPn/2). The Laplace transform jj,{s) = jliu = is) = 
/io(s/so)'^ is then derived by analytical continuation. Again, sq is an arbitrary 
reference frequency set to Sq = 1 s^^. As can be seen in Fig. |2l (Ar^(t)) can be 
approximated by a power-law in the range t G [0.015; 0.3] s. As a consequence, 
in the corresponding range s G [3; 60] s^^, the Laplace transform of the mean- 
squared displacement is well described by a power-law {Ar'^{s)) = a{s/so)~''- 
Fig. m shows the parameters a, b, fiQ and /5 as a function of aging time t^, aver- 
aged over a set of six different realizations, using the time rescaling described 
in section H?T1 

The exponent b decreases from about 1.9 ±0.05 to 1 ±0.05 while the exponent 
/3 increases from about zero at low aging times to 0.82 ±0.1 at the end of the 
experiment. This is consistent with the change from a viscous dominated to an 
elastic dominated behavior. The error bars on f3 take into account the depar- 
ture between the measured and extrapolated values, which are not systematic 
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(see another realization in (jl4j)). 



Bringing together our data for {Ar'^{s)) and fi{s), and using the modified 
Generahzed Einstein Relation (ITTl) . the GER ratio <9(s), parametrized by t^,, 
is given by : 

in the range s G [3, 60] s^^. 

After analytical continuation of Q{uj) = 0(s = ~iuj), we derive the effective 
temperature Teff(cij), at a given t^, from : 

W = (^M)'---^ (15) 

4fc/xo cos(p7r/2) 

in the range u G [3; 60] rad.s^^. In a system at equilibrium, we should have 
= T and 2 — 6-/5 = 0, which leads to 0(s) = T and Tofj(ct;) = T. This is 
not the case here for the aging colloidal glass. 

The dependence of the effective temperature T^s with aging time t^;, at dif- 
ferent frequencies, is shown in Fig. |5l The results have been averaged over six 
realizations. At the earliest ty^, the effective temperature is close to the bath 
temperature T = 300 K. Upon increase on t^, Tes increases up to 2-3 times 
the bath temperature and then decreases back upon further increase on tw 
Such a behavior - an increase of Tefj followed by a decrease - was observed 



for the first time in a colloidal glassy system in (Il4j ). where a scenario for 
this non-monotonic behavior was suggested. In the range u G [3; 60] rad. s~^, 
Tgfj increases with frequency for all aging times as seen in Fig. O In these 
experiments, the set-up does not allow to measure Tes at longer ty^. Indeed, 
beyond = 200 minutes, the mobility modulus becomes smaller than 10~^ m 
s^^pN^^. In this range, the optical tweezer is not powerful enough to induce 
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a detectable motion of the bead. 

Here we remain in the frame of power-law approximation. This is obviously 
not true in the full time range, but at least in the range [0.015; 0.3] s. Signif- 
icant deviations from a power-law are measured at large aging times and are 
currently under study. 



5 Discussion and conclusion 



As previously reported in (jl4j ). the dependence of T^g on can be interpreted 
through the change in the glass relaxation times with aging time. In colloidal 
glasses, dynamic light scattering and diffusive wave spectroscopy experiments 



provide the opportunity to 
loidal particles diffusion (129 



3robe the relaxation times associated to the col- 
; 1371 ). Fig. [H] shows a scheme of the resulting 
relaxation times distribution function -P(t), centred around two characteristic 
relaxation times Tfast and TsIow Upon increasing t^, part of the modes distri- 
bution function, centred around Tgiow, shifts to larger times, while the mode 
at Tfast ~ 0.1 ms remains unchanged. 

We suggest the following scenario for the dependence of the effective tempera- 
ture on the aging time. When probing the colloidal glass at a typical frequency 
/ = 1 Hz, three situations successively occur upon increasing tw At the earli- 
est tw, the glass relaxation times are small compared to the observation time 
scale Tfast.siow ^1// ^"^^ not play any role on this observation time. The 
glass is "at equilibrium" . In this experiment, the bead is small enough to be 
sensitive to molecular fluctuations. This ensures that it can thermalize with 
the bath when the glass relaxation times do not play any role on the observa- 
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tion time scale. In this case, the reading of the bead thermometer reduces to 
the bath temperature T. Upon increase on t^, relaxation modes of the order of 
the observation time scale Tgiow ~ V/ appear in the glass as seen in Fig. Ob. 
On this observation timescale 1//, the system is out of equilibrium : the mea- 
sured temperature becomes different from the bath temperature. Deviations 
from the Generalized Einstein Relation are thus observed. Finally, for very 
long tw, the first situation is recovered. The slow and fast relaxation processes 
of the glass do not play any role at the experimental observation timescale 
because Tgiow ^1// and Tfast ^1//. The bead thermometer, only sensitive 
to the fast relaxation processes at the molecular scale, again thermalizes with 
the bath. The measured temperature is then expected to reduce back to the 
bath temperature T. Finally, in this scenario, deviations are only observed in 
the regime where the observation time scale is of the order of a characteristic 



relaxation time of the glass, ut ~ 1 (jl4j ). These different situations are clearly 
identified in Fig. [5l 



As described in 



38l). the temperature of an object is measured by cou- 
pling it to a thermometer during a sufficiently long time interval such that all 
heat exchanges between thermometer and system take place and the whole 

n 

system equilibrates. In (!2i), the authors consider an oscillator of characteristic 
frequency Uq as a thermometer, weakly coupled to a glassy system in which a 
regime with small energy flows exist. The glass exhibits two time correlation 
scales, a fast one and a slow one, such that: Rsrit — t^,) = -^dt^C srii — t.^,) 
is satisfied at short times t — ^ tyj and RAcitw/t) = j^dt^CAcitw/t) is 
satisfied if tw/t = 0(1). Rst and Cst (or Rag and Cag) are respectively 
the stationary or aging parts of the response and correlation functions. T is 
the bath temperature and T* is the temperature associated to FDR when 
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t^/t = 0{l). 



If ujq is high enough so that the thermometer evaluates the fast relaxation and 
quasiequilibrium is achieved on the time scale I/cuq, its asymptotic internal 
energy density reads -Etherm = kT. The reading of the thermometer is the 
bath temperature T. If c^o is very low, the thermometer examines the system 
in its long time scales behavior and one finds -Etherm = kT* , where T* is 
identified as the temperature of the system on this observation time. Any 
small but macroscopic thermometer, weakly coupled to the system, is shown 
to play the same role as the above considered oscillator. The role of ujq is then 
played by the inverse of the typical response time of the thermometer. 

In our experiments, the bead plays the role of a thermometer which character- 
istic frequency ojq is set to the probing frequency u. Indeed, the bead inertia 
is negligible below cUmertia = 10^ rad.s^^ for the highest value of the fluid vis- 
cosity (~ 10 Pa.s). Experimentally, the observation time scale l/cu is much 
higher than the typical response time of the bead thermometer l/ct;inertia- As a 
result, the response time of the bead is set to the observation time scale 1/uj. 
It follows that by changing this observation time scale, the thermometer can 
explore or be sensitive to, different relaxation time scales of the glass (see Fig. 

ED. 

In models and simulations - where the characteristic relaxation time of the 
system is set to the waiting time - deviations are observed when the charac- 
teristic observation time 1/uj is of the same order or greater than the aging 
time tw Experimentally, deviations are observed in the other regime when 
l/uj <^ t^. In our sense, this discrepancy arises because the aging time 
is not the only relevant parameter to describe aging, and that microscopic 
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processes, characterized by the distribution of relaxation times in the system, 
must be considered. In our experiments, deviations are observed when ujt ~ 1. 
This suggests that, besides the waiting time t^, the distribution of relaxation 
times must be included in models to get an accurate description of aging. 



In the past, FDR have been experimentally investi gate d in the same colloidal 



13l ) . In the dielectric mea- 



glass with electrical and rheological measurements % 
surements, the effective temperature decreases with and tu, and reaches the 
bath temperature at high frequency and aging times. In the rheological mea- 
surements, no deviation from FDR could be detected. More recently, FDR, in 
the form of a Generalized Stokes-Einstein relation, has been investigated by 



combining dynamic light scattering and rheological measurements (|15i). The 
authors assume that the bulk stress relaxation probed by a rheometer is the 
same as the local stress relaxation affecting the probe particles. The effective 
temperature is shown to increase with frequency and aging time in the range 
tu, G [82; 135] min and uj G [20; 100] rad. s~^ for a 3% wt suspension. The 
authors suggest the existence of two regimes: ujt -C 1 leading to the bath 
temperature T and ujt ^ 1 leading to Teg, which is not consistent with our 
interpretation despite the authors claim. Other groups have recently investi- 
gated the deviations from FDR, by performing a test of the Stokes-Einstein 



relation, using optical tweezers (JIG 



171 ). In one case, they found that no de- 



viation could be detected in the aging colloidal glass over a wide range of 
frequencies / G [1.2; 12000] Hz. In the other one, deviations are observed at 
large aging times. Although the reason of these discrepancies have to be elu- 
cidated, our experimental findings are fully supported by our interpretation 
based on microscopic processes probed in light scattering experiments (see 
Fig. ED. 
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At this stage, whether these experimental investigations with the same observ- 
ables, in the same colloidal glass, are in contradiction still remains an open 



question ([l4 



15 



16l ). Furthermore, the definition of an effective temperature. 



with a thermodynamic meaning, implies that thi s q uantity is independent 



of the chosen observable, which is not verified in (jl3l ) where rheological and 
dielectric measurements do not lead to the same effective temperature. Is 
the effective temperature defined through a generalisation of FDR a relevant 
concept ? Does it depend on the chosen observable, as shown in the glass 



phase of Bouchaud's trap model (jld ). or not, as in simulations on a binary 
Lennard- Jones mixture jsol) ? More experimental investigations are needed to 



understand these contradictory findings (see (1401 ) for a review 



In conclusion, this work provides an experimental and theoretical investigation 
of the Generalized Einstein Relation in an aging colloidal glass. We interpret 
the observed deviations from GER as directly related to the change in the 
glass relaxation times. Deviations are observed in the regime ujt ~ 1. 



As a final comment, let us add that dynamical heterogeneities are recognized 
as a general feature of slow dynamics encountered in supercooled liquids and 



glasses (|41 
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43l ). In supercooled liquids, deviations from the Einstein re- 
lation were explained by such dynamical heterogeneities (141'). Up to now, in 
models and simulations, deviations from FDR have been found in homoge- 
neous glassy systems and interpreted in terms of an effective temperature. 
Moreover, the effective temperature is defined as an ensemble-averaged quan- 
tity in experiments, models and simulations. Understanding how the concept 
of effective temperature interplays with such heterogeneities is an open ques- 
tion. 
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Fig. 3. Modulus (a) and pliase (b) of tiie complex mobility fjL{u) = \fl{Lo)\e^'^^'^^ 
of the bead as a function of tw, for various frequencies of the applied force; (c) 
mobility modulus versus lo for various t^', from top to bottom, varies from 8 to 
191 min, for this realization. At low t^, is nearly independent of u. Upon 

increasing tyj, the modulus is well fitted by a power-law |A(^)| = Mo('^/'^o)^ with 
(jjQ = 1 rad.s^^ and /3 only depending on tyj. The full circles in (b) correspond to 
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Fig. 4. Top: Brownian motion parameters a and 6 as a function of aging time, where 
(Ar2(s)) = a{s/so)~'', with sq = 1 s~^; Bottom : Mobility parameters //q and /? as a 
function of aging time, where jl{oj) = ixo{u)/u)o)^ ^xp (— i/37r/2), with ojq = 1 rad.s"^. 
The results have been averaged over six realizations, and rescaled in aging time. 
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Fig. 5. Effective temperature T^s of the colloidal glass as a function of aging time tyj , 
measured at different frequencies. Upon increase on t^, Teff increases up to about 
2-3 times the bath temperature and then decreases back. The results have been 
averaged over six realizations rescaled in aging time. The typical errors bars on 
Tf.ff have only been shown on the data at / = 1 Hz for clarity, but all error bars 
roughly have the same magnitude. 
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Fig. 6. Scheme of the relaxation times distribution function P(r) in the glass of 
Laponite (typically 2.5% wt) at different aging times tw Upon increasing tw, part of 
the modes distribution, centred around TsIow, shifts towards larger times, while the 
mode at rpast ~ 0.1 ms remains unchanged. The arrow represents the measurements 
timescale 1/f. 
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